Transfer-matrix method

For the method in optics, see Transfer-matrix method (optics).

In physics and mathematics, the transfer-matrix method is a general technique for solving problems in statistical mechanics.

The basic idea is to write the partition function in the form

\mathcal{Z} = \mathbf{v}_{0} \cdot \left\{ \prod_{k=1}^{N} \mathbf{W}_{k} \right\} \cdot \mathbf{v}_{N%2B1}

where v0 and vN are vectors of dimension p and the pxp matrices Wk are the so-called transfer matrices. In some cases, particularly for cyclic systems, the partition function may be written more simply as

\mathcal{Z} = \mathrm{tr} \left\{ \prod_{k=1}^{N} \mathbf{W}_{k} \right\}

where "tr" denotes the matrix trace. In either case, the partition function may be solved exactly using eigenanalysis. If the matrices are all the same matrix W, the partition function may be approximated as the N power of the largest eigenvalue of W, since the trace is the sum of the eigenvalues and the eigenvalues of the product of two diagonal matrices equals the product of their individual eigenvalues.

The transfer-matrix method is used when the total system can be broken into a sequence of subsystems that interact only with adjacent subsystems. For example, a three-dimensional cubical lattice of spins in an Ising model can be decomposed into a sequence of two-dimensional planar lattices of spins that interact only adjacently. The dimension p of the pxp transfer matrix equals the number of states the subsystem may have; the transfer matrix itself Wk encodes the statistical weight associated with a particular state of subsystem k-1 being next to another state of subsystem k.

As an example of observables that can be calculated from this method, the probability of a particular state m occurring at position x is given by:

\mathrm{Pr_m(x)} = \frac{\mathrm{tr} \left[ \prod_{k=1}^{x} \mathbf{W}_{k} \mathbf{Pj} \prod_{k'=x%2B1}^{N} \mathbf{W}_{k'} \right]} { \mathrm{tr} \left[ \prod_{k=1}^{N} \mathbf{W}_{k} \right] }

Where Pj is the projection matrix for state m, having elements Pj_{\mu\nu} = \delta_{\mu\nu}\delta_{\mu m}

Transfer-matrix methods have been critical for many exact solutions of problems in statistical mechanics, including the Zimm-Bragg and Lifson-Roig models of the helix-coil transition, as well as the famous exact solution of the two-dimensional Ising model solution by Lars Onsager.

See also

References